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Abstract 

Wc consider the radiation of photons from quarks scattering on color- 
magnetic monopoles in the Quark-Gluon Plasma. We consider a temperature 
regime T ^ 2Tc, where monopoles can be considered as static, rare objects 
embedded into matter consisting mostly of the usual "electric" quasiparti- 
cles, quarks and gluons. The calculation is performed in the classical, non- 
relativistic approximation and results are compared to photon emission from 
Coulomb scattering of quarks, known to provide a significant contribution to 
the photon emission rates from QGP. The present study is a first step towards 
understanding whether this scattering process can give a sizeable contribution 
to dilepton production in heavy-ion collisions. Our results are encouraging: 
by comparing the magnitudes of the photon emission rate for the two pro- 
cesses, we find a dominance in the case of quark- monopole scattering. Our 
results display strong sensitivity to finite densities of quarks and monopoles. 

1 Introduction 

Creating and studying Quark-Gluon Plasma, the deconfined phase of QCD, in the 
laboratory has been the goal of experiments at CERN SPS and at the Relativistic 
Heavy Ion Collider (RHIC) facility in Brookhavcn National Laboratory, soon to be 
continued by the ALICE (and, to a smaller extent, by the two other collaborations) 
at the Large Hadron Collider (LHC). Dileptons and photons are a particularly in- 
teresting observable from heavy ion collisions, since electromagnetic probes do not 



interact with the medium after their production, thus carrying information about all 
stages of the evolution [T]. Discussion of dilepton production in heavy ion colhsion 
experiments at CERN SPS |3] and at RHIC [Ij and their comparison to theory can 
be found in [HI El [7], for a recent review see also e.g. [Hj and references therein. The 
main contributions to the production rates considered so far are hadronic decays 
in the mixed and hadronic phases of the collision (the so-called hadronic cocktail), 
and quark-antiquark annihilation in the QGP phase. After a long history of experi- 
mental studies of dileptons produced by charm decay, NA50/60 experiments finally 
concluded P that they do see QGP radiation, at intermediate dilepton mass range 
1-3 GeV, as predicted in [B |10]. 

However, the experimentally observed excess in dilepton production at small 
Pt and small invariant dilepton mass (below rUp ~ 0.77 GeV) remains a puzzle: 
the sum of all known contributions fails to explain the data by a large margin. 
Motivated by this, we search for additional, unexplored mechanisms which might 
contribute to the dilepton production rate. In particular, we will focus here on the 
role played by color-magnetic monopoles: we want to estimate the contribution to 
dilepton production from quarks which scatter on them. This methodological paper 
is our initial step towards an exploration of the subject: by no means we claim any 
resolution of the puzzle. 

The so called "magnetic scenario" for QGP [H] suggested that the near-Tc region 
is dominated by monopoles. More specifically, these authors suggested to look at the 
magnetic sector as a (magnetic) Coulomb plasma of monopoles in its liquid form. A 
line of lattice-based results has led to a very similar conclusion [12] . This scenario 
has met with initial success by providing an explanation of the low viscosity observed 
at RHIC [m Us] due to the large transport cross section induced by scattering on 
monopoles. 

Lattice monopoles are defined by the procedure [H] which locates the ends of 
singular Dirac strings by calculating the total magnetic flux through the boundary 
of elementary 3-d boxes. Since this depends on a certain gauge fixing, for decades 
sceptics kept the viewpoint that those objects are just unphysical UV gauge noise. 
Yet, many specific observables - e.g. monopole density - produced very reasonable 
and consistent results, apparently independent of the particular lattice parameters 
[121 [IS]- More recent results on monopole correlations [15j quantitatively support 
the Coulomb plasma picture of Ref. [llj, providing further reasons to think that 
monopoles are not artefacts but meaningful physical objects, present in the QGP 
as a source of a Coulomb-like magnetic field on which charged particles (quarks) 
can scatter. We work under the same assumption in the present paper. We are not 
advocating this magnetic scenario, but rather use it to estimate the contribution of 
radiation on monopoles. We will not need any assumption about monopole coupling, 
internal structure or correlations, only their density. 

Although we were motivated by the dilepton puzzle and will eventually aim at 
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solving it as a final goal, we start from the simpler problem of soft photon radiation 
during the collision process. The emission of real photons is a process which is closely 
related to dilepton production: the latter takes place through the emission of a 
virtual photon. In the QGP phase, the leading perturbative diagram is the Compton- 
like process {qg — > 57 and the crossing diagram qq — > •yg), while perturbatively 
subleading bremsstrahlung diagrams {qq —>■ qq'j) and LPM-type resummed effects 
are in fact equally important [16J. 

The classical trajectory of a particle with electric charg^ e in the field of an 
infinitely heavy monopole with magnetic charge g takes place on the surface of a 
cone. The static monopole approximation is valid for a regime of temperatures 
T ^ 2Tc, where they can be considered heavy, rare objects embedded into matter 
consisting mostly of the usual electric quasiparticles, quarks and gluons. The Lorentz 
force acting on the electric charge is proportional to the product of both couplings 
[eg). Thanks to the Dirac charge quantization condition, eg = 1 and thus it is not 
a small parameter and it is T-independent. 

As a first step towards the solution of the problem, in the present paper we 
compute the radiation from a non-relativistic electrically charged particle moving 
in the field of a monopole along the classical trajectory, ignoring back reaction. A 
full quantum and relativistic study is postponed for future investigation. Below we 
will discuss the applicability limits of our approximation. Therefore, our present 
calculation cannot address the actual phenomenological questions yet, but we can 
get an insight into how sizeable the effect of monopoles can be. To this purpose, we 
will compare our results with the parallel computation of photon emission rate in 
the process of Coulomb scattering of quarks, in the same approximations, regarding 
the Coulomb problem as a benchmark for comparison. 

Let us start with a "naive estimate" for the ratio of emission rates of quark- 
monopole vs Coulomb scattering of quarks: 

r^^i^gfv^f^uM 
111 ^ m Uq' 

First, the emission amplitude from monopoles is suppressed by the velocity v of the 
incoming electric particle, because the underlying scattering happens due to Lorentz 
rather than Coulomb force. Second, it obviously contains the density of monopoles 
hm ~ T^/ln^(r), which is smaller than the quark density ~ 

On the other hand, this rate is enhanced by the ratio of coupling constants. 
The numerator includes the product of the electric gauge couphng constant e and 
the magnetic one g: in the units we are using, the Dirac quantization condition 

-'^Here and below we call e the strong interaction coupling constant, using QED-like field normal- 
ization of the fields and 47r for consistency with the textbook material we use: note that = a^. 
The name g is reserved to color-magnetic coupling, while the electromagnetic coupling will be 
denoted as eem, again with = aem- 
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implies eg = 1 (actually h), while the electric Coulomb scattering is proportional 
to = ~ l/ln(T) -C 1. The In^(T) in the numerator to a significant extent 
compensates the smallness of the monopole density um ~ l/ln^(T). Although 
formally still decreasing at large T, as is the corresponding ratio of the contributions 
to viscosity [13], large angle and even backward scattering induced by monopoles 
may make this ratio numerically enhanced. The reduced mass fi = m/2 enters the 
Coulomb scattering problem whereas, in the limit of infinite monopole mass, the 
quark mass m appears in Newton law. There is an additional relative enhancement 
in favour of monopoles, which is due to the Casimirs of the qq and qq potentials. 
We will work out these factors in the following; they roughly bring a factor 1/4 
compared to the q — M scattering. For typical (thermal) velocity of 0.7, as = 0.8, 
and the ratio of densities 0.2 we end up with a relative effect of order 1/2. This 
is obviously a good start, suggesting it is worth to examine the problem in more 
details. 

We will see below that the above crude estimate is indeed correct when the veloc- 
ity of quarks is sufficiently close to one (we will still be using the non-relativistic ap- 
proximation); it is also correct in the ultrarelativistic limit. There is, however, a very 
important effect, which noticeably enhances the relative contribution of monopoles. 
The above estimate holds for a fixed impact parameter. For soft radiation, only emis- 
sion from large impact parameters is relevant (and this is the only region where our 
approximations are in fact valid). However, finite densities of quarks and monopoles 
in the plasma provide natural cutoffs for maximal impact parameters. Moreover, 
since the density of monopoles in our temperature regime is much smaller than 
the density of quarks, the relevant upper cutoff on the impact parameter for the 
Coulomb problem is much smaller and it leads to a significant suppression of the 
Coulomb-induced emission rate relatively to the rate due to scattering on monopoles. 

Our final results are very encouraging. We find that the soft photon emission rate 
from quark-monopole scattering could be as large as the mechanisms previously ac- 
counted for. Therefore, the problem calls for further and more detailed investigation 
at the full quantum level, which will be our next step in this project. 

The paper is organized as follows. In Section 2 we provide a brief overview 
of the classical motion of an electric charge in a Coulomb magnetic field: even if 
these results are well-known, a short summary is useful, since we will make large 
use of them in the following. In Section 3 we compute the photon radiation rate 
for an electric particle moving in the field of a static monopole. This section can 
be regarded as an extra chapter for Ref. [17], in which radiation rates for several 
trajectories were computed. In Section 4 we present an estimate of our effect in 
the case of quarks scattering on monopoles in the QGP: we use parameter estimates 
which are typical of the QGP produced at RHIC. We compare our results with those 
obtained in the Coulomb problem, and discuss the validity of our approximations. 
We draw our conclusions in the last Section (5), where we also indicate future 
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improvements. Appendix A supplements Section 3 by providing details of analytical 
computations for the emission rate. Appendix B presents a calculation of the quark 
density in the PNJL model, a result which is needed when we apply our calculations 
to QGP. 



2 Classical quark-monopole scattering 

We consider the classical, non-relativistic motion of a charge in an external field 
[iHl Uni EO] . A pointlike magnetic charge g is the source of a Coulomb-like magnetic 
field 

r 

B = 9-r (2) 
The equation of motion of an electrically charged particle e in such a field is 

m— = e.x5 = --xr; (3) 

the static monopole is located at the origin and the vector f defines the position of 
the electric charge (see Fig. [T]). In the following, we set c = 1 for simplicity. We 
also use the convention = a, and therefore eg = 1. 



In this process, the kinetic energy of the electric charge is a constant: 

E = — — = const., (4) 



as is the absolute value of the velocity vector v. There is no closed orbit in the 
charge-monopole system: the electric charge is falling down from infinitely far away 




egr 



Figure 1: The motion of an electric charge in the field of a magnetic monopole. 



onto the monopole, approaching a minimal distance b and being reflected back to 
infinity. This is evident from the trajectory 



(5) 



A special feature of such a motion is that the conserved angular momentum is differ- 
ent from the ordinary case: the absolute value of the ordinary angular momentum 
is conserved, but its direction is not constant. The generalized angular momentum 



J — [fx mv\ — eg- 



(6) 



is an integral of motion. The trajectory of the electric charge does not lie in the 
plane of scattering orthogonal to the angular momentum; the angle between the 

— * 

vectors J and f is a constant and the electric charge is moving on the surface of a 
cone whose axis is directed along — J with the cone angle 9 defined as: 



sin^ 



mvb 



cos 9 



eg 



(7) 



The velocity of the electrically charged particle is 



df 
di 



1 



J X f 



+ 



vH 



J X f 



{h/vty 



zr — X f+ Vrf (8) 



where the angular and radial components of the velocity vector are 

J V 



V ~ 2 ' 

mr 



Vr — 



(9) 



asymptotically we have 



t=±oo 



= 0, 



t=±oo 



(10) 



while at the turning point of the path we have 



t=o 



mb'^ 



0. 



t=o 



(11) 



The azimuthal angle ip as a function of time can be obtained by integrating the 
angular velocity: 

1 ''^t 

(f — arctan — . (12) 

sin6' 





y [GeV-'l 



y [GeV" 

(b) 

Figure 2: Example of electric charge trajectory in the field of a magnetic monopole. 
(a) X and z axes: the trajectories of the incoming particle (from t = —oo to t = 0) 
and the outgoing particle (from t = to t = +oo) overlap in this case, (b) y and 
z axes. Solid line: incoming particle (from t = — oo to t = 0). Dashed: outgoing 
particle (from t = to t = +oo). (c) Three dimensional trajectory. For all plots we 
use V = 0.5; b = 1 GeV-\ m = 0.3 GeV. 



An example of trajectory is shown in Fig. |2} 



The acceleration of the electric charge, a = Cj, follows from Eqs. (|3|) and (|8|): 



eg ^ ^ 
a = — -V X r 



eg 



v.r, X r) X r. 



(13) 



We recall that is directed along the z axis (the cone axis) so that: 

^ _ eg{v^). 



m r 

eg 



-irl + rl)] 



\/ {mvby + {eg)'^ sin 9 cos 9 [— cos ip, — sin ip, — tan 5] . (14) 



3 Radiation 

Following Ref. [TTj we define the intensity dl of radiation into the element of solid 
angle dQ as the amount of energy passing in unit time through the element df = 
R^dVL of the spherical surface with center at the origin and radius Rq. This quantity 
is equal to the energy fiux density {S = ^n) multiplied by df: 

tt2 

dl = —Rldn. (15) 
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The energy dSft^^, radiated into the element of sohd angle dQ in the form of waves 



with frequencies in the interval duj/2Ti is obtained from Eq. (15) by replacing the 
square of the field by the square modulus of its Fourier components and multiplying 
by 2: 

d£^ = ^-^Rldn^ (16) 



where 
and 



rin 



H^ = ikx (17) 

ikRo r+oo 



We work in the dipole approximation, namely we neglect retardation effects. 
This approximation is valid provided that v <^ c. In this case, the field can be 
considered as a plane wave, and therefore in determining the field it is sufficient to 
calculate the vector potential. The intensity of the dipole radiation is: 

rf/= %^sin2^rffi, (19) 

4:71 

where a is the acceleration of the electric charge. Replacing dQ = 2iTsin6d6 and 
integrating over 6 from to tt, we find the total radiation: 

/=^(ee™a)2. (20) 

The quantity dS^^ of energy radiated throughout the time of the collision in the 



form of waves with frequencies in the interval du/2Tr is obtained from Eq. (20) by 
replacing a by its Fourier component a^j and multiplying by 2: 

d£u = t; (eemdujf (21) 

In a standard planar collision, the total radiation dn^ in a given frequency interval 
du can be obtained by multiplying the radiation dSi_j from a single particle (with 
given impact parameter p) by the measure 2'Kpdp and integrating over p from pmin 
to pmax- In our case, the motion of the electric particle takes place on the surface 
of a cone, therefore we have to find the corresponding measure for our process. 

If we project the conical motion of the electric particle on a plane orthogonal to 
J, it is possible to define the standard impact parameter p for the planar motion 
[20] . We can define a position vector R: 



^^Jx{fx J) 



cos 9 cos 9 



f — J ( f ■ J 



(22) 
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which is the projection of r onto the plane perpendicular to J, times a factor 1/ cos 6* 
chosen so that R and r have the same length. At t = —oo we have: 



B 1 



cos 6* 



V — J i V ■ J 



(23) 



The mechanical angular momentum of the projected motion is the same as the 
conserved total angular momentum of the motion in the monopole field: 

J = mRxR; \J\ = mp\R_^\. (24) 

From the above Equation we get: 

V {mvbf + {egf 

P = 25 

mv 

from which it is easy to obtain that pdp = bdb: the measure over which we need 
to integrate dS^ turns out to be equal to bdb for the conical motion. Therefore, we 
can identify b as the real impact parameter for the process that we are considering 
{b is actually the minimal distance between the electric particle and the monopole, 
which is reached at t = 0): 

^''^ / 2Mb^. (26) 



du Ju . du 

The limits of integration for b are related to the specific scattering process we are 
dealing with. For example, if we consider a single quark scattering on a single 
monopole, we have b^ax — > C)0, while bmin can be identified with the size of the 
monopole core. In the problem we are considering, in which there is a finite density 
of monopoles, bmax turns out to be finite, as we will see in the next Section. 

We start by calculating the total radiation throughout the time of the collision 



following Eq. (20) 



X= / /(t)c?t = / \d{t)\'dt] (27) 



the acceleration components in the coordinate space are given in Eq. (14). We 
obtain: 

2 {eg)W _ 2 jegfy^b' 
^ ' 3 "'"m2r(t)6 3^"^m2(t;2t2 + 52)3 ^^^) 

so that: 

2 {egfy ^ dr _n {egfv 
~ 3^" mW (r2 + 1)3 - 4"^" m263 ' 
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the total radiation can be obtained in the following way: 

. = = ^\J^ ij^-T^)- (30) 

We now proceed by calculating the Fourier transform of the acceleration a in the 
case of quark-monopole scattering: 

_ {egf ^^ cxp [iOt] cos [^arctant] 

^ ^-^^ exp[i^.]sin[earctan.] 

mhi if + 1)= 



= / ..3/2 (33) 



where 



e = (^g)^ ^ = (34) 



For positive uj we obtain (see Appendix A): 



rU(e-l)U^U(e-l),-l,2u; 



^ -1 t^^!(p-^)!r(«i^-p + ^ + i) ^ ^ 



i exp (—a;) cos 



2 



-UJ] 



^A;!(p-A;)!r(^-p + fc + l) 

C+1 



X 



(36) 



X 2'=-2r I p - 
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(37) 



where p is the smallest integer number larger than + 3)/2. p — 2 is the number 
of full rotations around z axis. U {a, b, z) is the confluent hypergeometric function 
with integral representation: 



U{a, b, z) 



Via) 



dt 



(38) 



and K„{z) is the modified Bessel function of the second kind: 



cost dt 



'o (t2 + 22)"+l/2- 



(39) 



Equations (35)-(36) are strictly valid for u; > and for any value of ^, except odd, 



integer numbers. When ^ is an odd integer number, the above formulas vanish 



identically, and the integral is given by Eq. (65) (see Appendix A). The behavior of 
{ax)u}^ {ciy)uj and (a^)^^ as functions of u and b is shown in the two panels of Fig. [3| 
The component {ay)^^ is purely imaginary and vanishes at = 0. The values of the 
acceleration components at = are: 



2v 
1 



cos 



[a. 



y)uj=0 



0: 



mb^ 



(40) 



The subleading small uo asymptotic behavior can be found in Appendix A. The 
photon emission rate is finite as — 0. It is of course how it should be: the 
corresponding number of photons dN^ = dl^/hu would show standard logarithmic 
IR divergence. 



4 Application to QGP 

In this Section we give a rough estimate of the effect that we are describing, in 
the case in which the electric charge is a quark q (or an antiquark q), scattering 
on a color-magnetic monopole in a deconfined medium. The medium contains a 
finite density of quarks and monopoles. We will compare our result to the radiation 
produced from Coulomb scattering of qq, qq and qq pairs. We consider a regime of 
temperatures ^ 2 Tc where, according to the magnetic scenario proposed in [TT] . 
monopoles can be considered as heavy, static particles. 

A quark moving in a deconfined medium acquires a thermal mass due to its 
interaction with the other particles of the medium. Lattice results for this quantity 
are available for three values of the temperature ^11 [22]. At T = 1.5Tc and T = 3Tc 
they find m = 0.8T, while at T = 1.25Tc they obtain m = 0.77T. At T = 2T^ 
we therefore assume a value of m ~ 0.8T, namely m ~ 0.3 GeV. For temperatures 
T ~ 2Tc, we can assume that quarks move with an average velocity w ~ 0.5 — 0.7. 
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/ (a.vOw 




/ Im[(flv)„]. 
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cj [GeV] 
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Figure 3: (a): Fourier components of the acceleration as functions of uj, for m = 0.3 
GeV, 6 = 1 GeV~^, v = 0.7. (b): Fourier components of the acceleration as functions 
of b, for m = 0.3 GeV, cu = 0.1 GeV, v = 0.7. 
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Figure 5: A charge scattering on a 2-dimensional array of correlated monopoles 
(open points) and antimonopoles (closed points). The dotted circle indicates a 
region of impact parameters for which scattering on a single monopole is a reasonable 
approximation. 



In Figure |4] we show the integrand for the total radiation in a given frequency 
interval, namely 27TbdSuj/du, as a function of uj and b. We need to integrate this 
quantity over the impact parameter. There is no upper limit on b in the case of 
scattering of one quark on one monopole. However, in matter there is a finite 
density of monopoles, so the scattering issue should be reconsidered. A sketch of 
the setting, assuming strong correlation of monopoles into a crystal-like structure, 
is shown in Fig. |5} A "sphere of influence of one monopole" (the dotted circle) gives 
the maximal impact parameter to be used 



The same is true for the quark-quark and quark-antiquark Coulomb scattering 
to which we will compare our results. The monopole density as a function of the 
temperature for SU (2) gauge theory has been evaluated on the lattice [15] . In order 
to account for the transition from SU{2) to SU{?>) gauge group, we scale these 
results by a factor 2: in SU{2) there is in fact one monopole species while in SU{3) 
there are two, identified by two different U{1) subgroups. The monopole density at 
T ~ 2Tc is um — 0.02 GeV^, which gives bmax — 1.8 GeV~^ The lattice also gives 
us information about the monopole size [23j: they turn out to be very small objects, 
having a radius tm — 0.15 fm=0.78 GeV~^. This is the bmin that we will use in our 
integration. 

We therefore obtain, for the total energy radiated in unit volume throughout the 
time of the collision in a given frequency interval: 

dE 2dK^ 2 2 ^ /'^™"%,^ ,2 „ 

— = --j—ngUM = -n^nM—OLem^ir / b\a^\ db (42) 

duj 9 du 9 671 h - 
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where the factor | = + | + comes from the different electric charges fur u, d 
and s quarks (the density of quarks is the sum of the densities of u, d and s 
quarks). 



4.1 Comparison with Coulomb scattering 

In this Section we give an estimate of the radiation produced in the scattering of 
gg, qq and qq pairs in the plasma. In the case of an attractive interaction between 
particles (namely in the singlet channel for the qq scattering and the antitriplet 
channel for the qq and qq scatterings), the formula for dE^jduj reads |17j : 



duj 
where: 



V 



y^em ) 1 

mi 



ni2 



12 e^-l 

+ TT- 



Q ■ 



1 + 



mira2 
mi + m2 



and H^}}(ii'e) is the Hankel function of the first kind: 



(43) 



(44) 



JJz) +iYJz). 



(45) 



In Eq. (43), {eem)i, and (eem)2, ""^2 are the electric charge and mass of the two 
colliding particles, while a = C^ag is the strong coupling constant multiplied by 
the corresponding Casimir factor for the channel under study: 

C^ = 4/3; = 1/6; = 2/3 ; C^ = l/3. 

When the interaction is repulsive (namely in the octet channel for the qq scattering. 



and in the sextet channel for the qq and qq scatterings), Eq. (43) gets modified as 
follows: 



dSu, 
doj 



l-Ka^Up- ( (eem)l {eem)2 



+ 



=2-1 



mi 



1712 



exp [— 27rz/] 



(46) 



We consider quark matter with three equal mass light flavors, namely mi = m2 = 
m; fi = 171/2. The total density of quarks and antiquarks can be obtained for 
example from the PNJL model (see Appendix B). 



Uq ~ 2.8T^ ~ 0.12 GeVl 



14 



The total energy radiated in unit volume throughout the time of the collision in the 
case of Coulomb scattering can be obtained from the following formula: 



3f4 



-a. 



f 



4 



-a. 



-a. 



f 



-a. 



exp[— 27rz/as/6] 



+ 3 I -a, 



f 



+ 6 



where 



/(«)= H^^'(^ue 



2 

+ - 



exp[—27ii'as/3] 



bdb. 



1 



(47) 



(4J 



Eq. (47) is the total radiated energy, it takes into account all possible color chan- 
nels for qq, qq and qq scatterings, and all possible flavor combinations, bmax can 
be estimated directly from the quark density: in our temperature regime, it turns 
out that bmax — 1 GeV^^. Our results for the ratio of the total energy radiated 
throughout the collision time in the case of quark-quark and quark-monopole scat- 
tering are shown in Figs. [6j The left panels show this ratio for bmax ~^ oo, while 
in the right panels bmax is finite and fixed by the corresponding densities. This is 
useful to understand how strongly a finite bmax influences our results: it turns out 
that the cutoff dependence of dll/dw is dramatic. Without cutoff, this quantity is 
much larger in the case of a Coulomb scattering, while the opposite is true in the 
case of a finite bmax- This effect is qualitatively true for all values of v that we have 
considered, but obviously the relative magnitude for qq and qM scatterings depends 
on the specific value of v that we choose. 



4.2 Validity of our approximations 

Our results are obtained through a series of approximations: 

• non-relativistic approximation: this is obviously valid if the velocity v of quarks 
is not too large compared to the speed of light, namely f -C 1; 

• dipole approximation vs retarded emission. It is valid if the radiation wave- 
length is large. The retardation effects can be neglected in cases where the 
distribution of charge changes little during the time a/c, where a is the order 
of magnitude of the dimensions of the system. This is true if f <^ 1, which 
coincides with the condition for the non-relativistic approximation; 

• classical trajectory is used, without back reaction of radiation: in all our for- 
mulas, we assume that the particle moves along a trajectory which is the 
solution of the classical equations of motion. This means that the energy lost 
by the particle through the radiation process is negligible. This approximation 
is valid in the limit a; <C m/e^; 
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Figure 6: Left column: ratio of dE/dw for Coulomb and quark-monopole scattering. 
In both cases, the integral over b is taken up to oo. For these plots we use ag — 
0.8, m = 0.3 GeV and: (a) v = 0.3, (c) v = 0.5, (e) v = 0.7. Right column: same 
as in the left column, but the integral over b is taken up to the corresponding bmax- 
For these plots we use as = 0.8, m = 0.3 GeV and: (b) v = 0.3, (d) v = 0.5, (f) 
V = 0.7. 
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• classical approximation: naively, the emission of soft photons can be described 
classically for 



Within the same approximation we can ignore the recoil effect (energy-momentum 
conservation). 

The full quantum treatment of the radiation should include the back reaction 
of the radiation, one has to evaluate the nondiagonal matrix element of the 
dipole moment between the initial and final scattering states, with different 
energies. Such quantum states for quark-monopole problem were found in 
|24j and recently for gluon-monopole problem in [13]. However, the matrix 
elements have not been computed yet. Since it was done for quantum Coulomb 
scattering, by A.Sommerfield in 1931, we can use those results in order to have 
at least some qualitative estimate for the accuracy of classical description. 

We compare the total radiation dn^) in a given frequency interval duo, in the 
case of scattering on a single particle, namely we integrate over the impact 
parameter h from to oo. In the classical case we have: 



Hl'/M (49) 



did 3m^v^ 
which is to be compared to the quantum expression [25j : 
dK{uj) _ 2 647r2 p' 1 1 ( d 2 



du """"^ 3 p ij) - p'Y (1 - e-2-/3') (e2-/3 - 1) V di 

(50) 

where 

13 = — = — mv =2p mv' =2p' p' = \/ ~ ""^^ 



and 



F{0 = 2F,itl3',t(3; 1;0 ^ = - 



4pp' 



2-Fi(a, b; c; d) is the Hypergeometric function. The two curves corresponding to 
classical and quantum scatterings are shown in Fig. [7} It is thus evident that 
within the energy region plotted the classical result is a very good approxima- 
tion of the quantum one. Note that there is a maximum u, which is equal to 
the energy of the incoming particle beyond which due to energy conservation 
the quantum formula are not applicable. 
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Figure 7: Comparison between classic and quantum result for dK^/duj as a function 
of uj, with as = 0.8, v = 0.7, m = 0.3 GeV. 



• two body vs many body scattering. So far, we have limited our analysis to two- 
body scattering, while mimicking the effects of multiple interactions by finite 
quark and monopole densities and the maximal impact parameter. The strong 
dependence on the maximal impact parameter cutoff observed in our results, 
reflects additional shortcomings of our approach. The origin of this problem is 
obvious: on the one hand soft radiation is emitted from large distances; on the 
other hand, too large distances are precisely governed by multiple scattering. 



5 Conclusions and outlook 

The present paper is the first step towards understanding whether the contribution 
of the quark-monopole scattering process in QGP is or is not important for photon 
and dilepton production. Our purpose at present was pretty modest: we wanted 
to evaluate the magnitude of soft photon radiation from this process and compare 
it to the one produced in Coulomb quark-quark scattering. A qualitative estimate 
outlined in the Introduction, said that while the monopole density is small at high T, 
nm ~ (1/ ln(T))^, the square of in the electric scattering cross section compensates 
two out of three such logarithms. So, parametrically, the process we consider is 
subleading at very large ln(T). On the other hand, the charge-monopole scattering 
at large angles and small impact parameter tends to be much larger than the charge- 
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charge one: we found that this enhances radiation, similarly to how it worked for 
transport processes fl^. This represents an encouraging starting point for future 
work. 

We have calculated the photon radiation rate for quarks scattering on monopoles 
in a thermal medium which contains a finite density of both particles. We worked 
in the classic, non-relativistic approximation and neglected retardation effects and 
back reaction. Therefore, our calculation has a rather methodological status, it 
cannot address the actual phenomenological questions, such as the experimentally 
observed excess in dilepton production at small pt and invariant mass below rup. We 
need to improve the present paper in many directions: first of all, a full quantum 
and relativistic calculation will be performed, also taking into account back reaction 
of the radiation. This can in principle be done using non-diagonal matrix elements, 
calculated between quantum scattering states such as those which were found in 

[iniEi]. 

Then we need to take into account the fireball evolution, in order to be able to 
quantitatively compare our results to the experimental data. Dramatic expansion 
of the fireball and long duration of the near-T^ phase leads to the conclusion that 
soft dileptons currently constituting the puzzle come from the end of the evolution 
of QGP, T ~ ITc, rather than the beginning of it, in the temperature regime we 
discussed above. Unfortunately, the near- Tc region is very complicated and quite 
challenging theoretically. In this region monopoles become lighter and thus dynam- 
ical and relativistic, matter is no longer an electric near-perturbative plasma with 
at least well-defined counting rules, but a strongly coupled liquid made of all kind 
of quasiparticles. We hope to address those issues elsewhere in our future work. 
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Appendix A 

In this Appendix we explicitly calculate the Fourier transform of the a components. 
We start from (a^): 

/ N _ (e^)^ n [^^^] [^arctant] 
mWvi (t2 + 1)3/2 

If we set t = i + ir, dt = idr. Both the square root and the arctan have branch cut 
singularities from r = to r = oo and from r = — oo to r = —2. This is clear if we 
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write the relationship between arctan and log: 

arctan {i + ir) = iarctanh (1 + r) = i log ^ — — ^ . (52) 

We calculate the integral along the contour C shown in Fig. [8j The integral can be 
split into four different contributions: 

• the small circle 7 of radius e 

• the upper and lower segments around the branch cut 

• the big circle F 

We have 

exp [iut] cos [^arctant] 



(t2 + l 



,3/2 



-dt 



r n f /"°° /" f f^ exp [-lut] COS [^aictan (i + ir)] , 

^,expi-^i|y+y^+yjj (U.-.^)^/- — '''' 

Along the upper segment, where r — > r + ze, we have: 



2^3/2 _ , , „2^3/2 

ZTT (54) 



(_2r-r2)'/' = -^(2r + r2) 

2 + r 



log ( -^-^ 1 = log 
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while along the lower one we have 



2\3/2 



log 



2 + r 



= log 



2 + r 



r 



ITT. 



(55) 



Therefore we can write: 

exp [—Cot] cos [^arctan (i + ir)] 



OO /■£ 
+ 



(_2r - T^f' 
exp [-cDt] cos (log I ^ I + «7r)] 



(56) 



i (2t + t2)='/' 



(ir + 



exp [—cot] cos (log 



I 2+r I 

I T I 



in)] 



i (2t + t2)='/' 



dT 



2i cos 



2 



exp [ — '^'T\ COS ] ^ iO^ 



(2t + t2)^/2 



dT 



where we used the relationship 



a + /5 

cos q; + cos p — 2 cos — - — cos — - — . 



We therefore have: 

exp [— i"r] cos [(arctaii (/ + ir)] 



OO fe 
+ 



(-2t-t2) 



2^,3/2 



dT 



I COS 



2 



exp [—cur] 



2 + r 



-m 



+ 



(2r + r2)='/2 

The above integral contains two terms. The first one gives (in the limit e ^ 0): 



(57) 
(58) 

dT. 



I COS 



2 



r exph 
io (2t + T 



-o'r] /2 + r 

2)3/2 



= I COS 



e7r\ 1 



-?/2 
dT 



(e-1) 



(e-i),-i,2^ . 



(59) 



The second term needs to be integrated by parts p times, where p is the smallest 
integer number > All boundary terms are either divergent or vanishing as 

e — > 0. The divergent ones are exactly cancelled by corresponding divergent terms 
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coming from the integral over the small circle 7. The finite part is: 



exp(-a;r)(2 + r)^^-')/' _ 4r - p) p\ ^ (-c^)^ 

-dr - 2^ . . X 



X / exp(-(I;T)[2 + T](^-=^)/'-(^-'=)Tf-(^+^)/2dr= (60) 
Jo 

_4r(5±!!_p)p! P 



E 



x2''-^r {p-^]u (p-^,k- l,2u 



The above contribution vanishes for odd, integer values of ^. The contribution 
coming from F vanishes identically, while we get a nonvanishing contribution from 
the integral over the small circle 7. This contribution is divergent for noninteger 
or even ^, exactly cancelling the divergent contribution coming from the segments 
along the branch cut. When ^ is an odd integer number we get a finite contribution. 
We redefine t = i + ieexp{i9), dt = —eexp{i9)d9 and get: 



exp [—cut] cos [^arctan (i + ir)] 

or 



-2r - T 



2^3/2 



^^-0, f^^e'^^"-^ cos [^arctan (i + ie exp(i^))] 



L 

1 (-2eexp(i^) - e2 exp(2i^))^/^ 



(-2eexp(ie) - e2exp(2i^))^/^ 

27r „ie„-e£i)e'^ r ^ / ^gi 

:log 



— ee / ^ cos 



2i ° V 2 + ee*' 



dQ 



2- ,-'^/2,, 



2(-e)V2yo (2 + eexp(i^))^/^ 



2 + ee*^ / ^ V 2 + ee 



dQ 
(61) 



For e — > 0, the first term in the parenthesis gives a finite contribution only for ^ = 1, 
a value which is never reached in practical cases, as we will see. For ^ > 1 its 
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contribution vanishes identically. The second term gives 



uj p2'K ^—i9/2^—eu)e 



ee 



de 



2. g-2^eg-e^e'« (^2e~'' + 6^^'^^^^ 



r)(5 + l)/2 



-de 



2n 



2(_e)(5+i)/2 



,n 2^e ^ ^"(-e)"e-^ ^^^^^ (2e-^)'=e«-3)/2-'^(i^)! 
2^ - 



n=0 



fc=0 



fc=0 



(A; + 2)! k\{^-k)\ 



' ^,+2 (_l)/c-(€-3)/22/c(i^), 

-27r 2^ 



fc=0 



(A; + 2)! k\{^-^-ky. 



(62) 



where we have taken into account the fact that the finite contribution from the sum 
over n comes from n = k + 2. The final result for (a^)^ is therefore 



+ 



^2Jj2y^ 

k 



exp {—00} cos 



2 



1„ fl 



4 V 2 



T[-{i-l)]U{-{i-l),-l,2Q 



-uo] 



^A;! {p-k)\V{^-p + k + l) 

e + 1 



(63) 



X 2'=-2r I p 



U \p- ^^,k~l,2uj 



(64) 



for any value of ^ other than odd-integer, and 



(ctx). 



(e(?)2 e^^_^^^/' uj^+^ (_i)fc-{€-3)/22fc(i_3). 



m?lP'V^ 2 



-2. ^ 



A:=0 



(A; + 2)! A;!(S^-A;)! 



(65) 



for odd-integer ^. 

The component (ay)t^ slightly differs from {ax)u)'- the sin integration gives a purely 
imaginary contribution. The second term in Eq. (63) has a minus sign. 
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The first terms of tlie components expansion around 00 = liave tlie following 
assymptotic behavior: 



2v 

[axju^^o ^ y cos 



+ UJ : COS 



+ 



-Harmonicnumber ( p 



e + 3 



Harmonicnumber 



e + 1 



(66) 



+ \n4:-2{^'^ -l)\n 



-UJ 



^2 _ I 



k=3 



k\{p-ky.T [^-p + k + l] 



{ay)^^o ~ 2«6cos 



UJ 



{eg)v 



a/ {egY + {mvby 



UJ 



b^eg (-1 + 27 + 2 In + 2 In {^)) 
Iv^iegY' + {mvby 



where Harmonicnumber (2;) = + 1) + 7 and 7 is the Euler's constant and "^{z) 
is the logarithmic derivative of F-function. 



Appendix B 

In this appendix we briefly recall some aspects of the PNJL model [26], which we will 
then use to calculate the density of quarks in the medium. This model successfully 
describes QCD thermodynamics in the temperature regime we are interested in, by 
coupling quarks to the chiral condensate and to a temporal background gauge field 
related to the Polyakov loop. The Euclidean action of the three-flavor PNJL model 
is 

Se{^P, V^ = J^^ ^'^ dr d^x [^t d^^ + u^ilj, ilj\ 0)] - ^ U{^, T). (67) 

Here Ti is the fermionic Hamiltonian density given by: 

n = -iip^ (a ■ V + 74 mo - 0) ^ + V(?A, ^^) , (68) 

where ip is the Nf = 3 quark field, a = 7o 7 and 74 = 270 in terms of the standard 
Dirac 7 matrices and mo = diag{mQu,mod,mos) is the current quark mass matrix. 
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V('0,'0^) contains two parts: a four-fcrmion interaction acting in the pseudoscalar- 
isovector/scalar-isoscalar quark- ant iquark channel, and a six-fermion interaction 
which breaks Ua{1) symmetry explicitly: 



+ 



G 
K 



f=u,d,s 



det {ipi (1 + 75) ipj) + det {ifji (1 - 75) ipj) 



(69) 



Quarks move in a background color gauge field 4> = A4 — iAo, where = 
Sf.ogA^t'' with the SU{3)c gauge fields and the generators = A"/2. The 
matrix valued, constant field (j) relates to the (traced) Polyakov loop as follows: 



-Tr 



Pexp i / drA^ 



-Tre*<^/^. 



(70) 



The thermodynamic potential of the system is: 



(T, /.) = W ($, T) + ^ + 1^ - ^<,a, 

13 



(71) 



2^ ^^'^3 |ln "1+ $e-(^^-^-^^)/^ + $e~^K^-''^)/^ + e~^K^-''^)/^] 



where 

/ K 

2G{ilJii>i), nj = -nj^ Epj = ^p'^ + mj, m, = moi-ai-—ajak. 

(72) 

By minimizing the thermodynamic potential one can obtain the behavior of the 
chiral condensates o"j and of the Polyakov loop $ as functions of the temperature 
and chemical potential. After this procedure, one is able to evaluate many thermo- 
dynamic quantities. For example, the quark density we are interested in is the sum 
of the densities of quarks and antiquarks and can be obtained through the following 
formula: 

dQ dQ 

= + ^ ^3 
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Figure 9: Total quark density as a function of the temperature (PNJL model result). 



and it has the following explicit form: 



Nf 



n„ 



3 exp[/i/T] (exp[2/i/T] + $ exp[2Ep/T] + 2$ exp[{Ep + fi)/T]) 
_exp[3Ep/T] + exp[3/i/T] + 3$ exp[(2Ep + fi)/T] + 3$ exp[(Ep + 2/i)/T] 



3 (1 + <^exp[2{Ep + /i)/r] + 2$exp[(igp + /i)/r]) 1 
l + exp[3(Ep + /i)/T] + 3$exp[2(Ep + /i)/T] + 3<l>exp[(Ep + /i)/T]J ^ ^ ^ ^ 

We show the behavior of quark density in Fig. [9} for typical PNJL model parameters 
taken from the last paper of Ref. [26] . 
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